ApJ accepted [ 1 8 April 2007] 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



THE DISTRIBUTION OF ACTIVE GALACTIC NUCLEI IN CLUSTERS OF GALAXIES 

Paul Martini 

Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, martini@astronomy.ohio-state.edu 



o 
o 

(N 

a 

>—> 
o 

6 



(N 
> 

in 
in 

m 

o 
l> 
o 



X 



John S. Mulchaey, Daniel D. Kelson 

Carnegie Observatories, 813 Santa Barbara St., Pasadena, CA 91 101-1292 
ApJ accepted [18 April 2007] 

ABSTRACT 

We present a study of the distribution of AGN in clusters of galaxies with a uniformly selected, spectroscop- 
ically complete sample of 35 AGN in eight clusters of galaxies at z = 0.06 — > 0.31. We find that the 12 AGN 
with L x > 10 42 erg s" 1 in cluster members more luminous than a rest-frame M R < -20 mag are more centrally 
concentrated than typical cluster galaxies of this luminosity, although these AGN have comparable velocity 
and substructure distributions to other cluster members. In contrast, a larger sample of 30 cluster AGN with 
Lx > 10 41 erg s" 1 do not show evidence for greater central concentration than inactive cluster members, nor 
evidence for a different kinematic or substructure distribution. As we do see clear differences in the spatial 
and kinematic distributions of the blue Butcher-Oemler and red cluster galaxy populations, any difference in 
the AGN and inactive galaxy population must be less distinct than that between these two pairs of popula- 
tions. Comparison of the AGN fraction selected via X-ray emission in this study to similarly-selected AGN in 
the field indicates that the AGN fraction is not significantly lower in clusters, contrary to AGN identified via 
visible-wavelength emission lines, but similar to the approximately constant radio-selected AGN fraction in 
clusters and the field. We also find significant evidence for variation in the AGN fraction between clusters and 
explore the dependence of cluster AGN fraction on redshift, velocity dispersion, amount of cluster substruc- 
ture, and fraction of Butcher-Oemler galaxies. While we see weak evidence for several trends, there are several 
correlations between these four parameters in our small sample of eight clusters that preclude identification of 
which one(s) most strongly influence the cluster AGN fraction. 

Subject headings: galaxies: active - galaxies: clusters: general - galaxies: evolution - X-rays: galaxies - 
X-rays: galaxies: clusters - X-rays: general 



1. INTRODUCTION 

Early work on clusters of galaxies found that emission-line 
galaxies and active galactic nuclei (AGN) were rar er in cluster 
members than in regions of lower galaxy density (lOsterbrockl 
Il960t lGisler|[l978t iDressler et al.lll985l) . Extensive work on 
cluster galaxy populations demonstrated that cluster galaxies 
were dominated by early-type galaxies with old stellar popu- 
lations and the fraction of the se quiescent ga laxies increased 
with galaxy surface density (Dressier 1980). These differ- 
ences are now ascribed to multiple physical mechanisms that 
cause cluster members to possess less of the cold gas neces- 
sary for young star formation and accretion onto supermassive 
black holes. 

The relative rarity of AGN in clusters of galaxies made it 
difficult to acquire suffi ciently large samples f or demographic 
studies. For example, IDressler et al.l ([1985) obtained spec- 
tra of over 1000 cluster members in 10 low -redshift clusters 
of galaxies and identified AGN in only 1% of them. This 
low fraction discouraged large spectroscopic studies whose 
primary purpose was to identify AGN in clusters of galax- 
ies. However, the advent of the Chandra X-Ray Observa- 
tory has provided a new means to efficiently identify AGN 
in clusters. Motivated by early evidence of an excess of X- 
ray point sources in the fields of rich clusters of galaxie s 
(Cappi et al. 200lt jSun & Murrayll2002T: iMolnar et alJl2002l) . 
we (iMartini et al.ll2002l) obtained spectra of the bright coun- 
terparts to X-ray sources in the field of the rich cluster 
Abell 2104 (z = 0.154). This study identified six bright cluster 
galaxies coincident with luminous X-ray sources (Lx > 10 41 



erg s" 1 ), or approximately 5% of cluster members more lumi- 
nous than Mr = -20 mag. While only the most X-ray lumi- 
nous of these six galaxies have the emission-line diagnostics 
of AGN, the remaining X-ray sources are more likely pow- 
ered by accretion onto a supermassive black hole than other 
plausible sources of lower-luminosity X-ray emission, such as 
a population of low-mass X-ray binaries (LMXBs), hot halos 
of diffuse gas, or star formation. 

Motivated by the high AGN fraction in Abell 2104, we 
extended our survey to seven additional, low-redshift clus- 
ters of galaxies (z = 0.06 — > 0.31) and found that approxi- 
mately 5% of galaxies more luminous than Mr = -20 mag 
host AGN more luminous th an Lx = 10 41 erg s " 1 in the 
broad (0.5-8 keV) X-ray band (iMartini et al.ll2006l hereafter 
Paper I). As was the case for Abell 2104, we found that 
most of the cluster galaxies with X-ray counterparts did not 
show obvious AGN spectral signatures in visible-wavelength 
spectra. As LMXBs or thermal emission from hot halos 
have also been observed to produce luminous X-ray emis- 
sion, particularly from bright, early-type galaxies, we used 
the multiwavelength spectral shape of these sources to deter- 
mine if they were AGN. Comparison with relations between 
X-ray and B-band lumi nosity for local early- type galaxies 
dominated by LMXBs ((Kim & Fabbiano 2004) and hot gas 
dO' Sullivan et all 120031: ISun et al.l 12005b showed that these 
AGN candidates were on order 1-3 orders of magnitude 
more X-ray luminous than expected from relations based on 
these other sources of X-ray emission. We therefore con- 
cluded that most of these sources were AGN. Only a small 
fraction remain consistent with other sources of X-ray emis- 



2 



MARTINI, MULCHAEY, & KELSON 



sion and they are not included in the present study. 

The goal of seeking a larger sample of AGN in clusters of 
galaxies was to use the AGN population to explore the mech- 
anisms responsible for fueling AGN, motivated by the simi- 
lar use of the cluster environment to explore other aspects of 
galaxy evolution. Comprehensive studies of clusters of galax- 
ies have shown that the cluster environment contains a pro- 
foundly different distribution of galaxy populations from the 
field. The morphology-density relation expresses the obser- 
vation that lenticular and then elliptical galax ies are dominan t 
in progressively richer galaxy environments (Dressier 1980). 
The amount of current star for mation in galaxies also declines 
toward the center of clusters (Fish eret alJl!998l ). Galaxies at 
the centers of rich clusters tend to be dominated by old stellar 
populations and have no active star formation, while a pro- 
gressively higher fraction of poststarburst and starburst galax- 
ies are found toward the outskirts. 

Kinematic studies of cluster galaxies show that the popu- 
lation of galaxies with current or recent star formation have 
the highest velocity dispersions, poststarburst galaxies are 
intermediate, while the galaxy population dominated by a 
passive stellar popul ation has the lowest velocity dispersion 
(Dressier et al. 1999). The higher velocity dispersions (and 
greater radial extent) suggest that galaxies with current or re- 
cent star formation are less virialized than more passive galax- 
ies. These galaxies with young stellar populations may on av- 
erage have entered the cluster more recently and have higher 
line-of-sight velocity dispersions because they remain on pri- 
marily radial orbits. 

The distribution of AGN in clusters of galaxies could pro- 
vide similar information on the origin of AGN in clusters. In 
particular, their distribution is a valuable test of the standard 
paradigm for AGN fueling, namely the m erger of two gas-rich 
galaxies ("e.g jBarnes & Hernquistll 992). The low AGN frac- 
tion in clusters is commonly ascribed to the lower merger rate 
in clusters due to the high velocity dispersion that precludes 
the formation of bound pairs, in spite of the high galaxy den- 
sity, and the lower fr action of galaxies with subst antial reser- 
voirs of cold gas (e.g. lGiovanelli & Havnesf l985). If this pic- 
ture is correct, then AGN should be more common at the out- 
skirts of clusters where members are relatively rich in cold 
gas, as well as in lower velocity dispersion clusters. Many 
galaxies also enter the cluster potential in low velocity disper- 
sion groups that may produce a relative increase in the AGN 
fraction at larger distances from the cluster center or associ- 
ated with distinct substructure within the clusters. These rela- 
tively recent entrants into the cluster potential would also have 
not yet virialized and could have a larger velocity dispersion 
than the old, passively-evolving galaxies at the center of the 
cluster potential. 

In s tudies of lower velocity dispersion groups, She n et all 
d2007l) found AGN in ~ 7% of galaxies at z ~ 0.06, yet these 
AGN were only identified in visible-wavelength spectroscopy 
and were not detected in their X-ray (XMM/Newton) observa- 
tions and must have Lx < 10 41 erg s . This result suggests 
that the X-ray and visible-wavelength properties of typical 
AGN in lower-density environments may be different from 
the typical AGN in higher-density environments. Previously, 
iBest et alj (|2005|) showed that the fraction of AGN selected 
by emission-lines from SDSS slightly decreases for galaxies 
with a larger number of luminous neighbors, yet the fraction 
of radio-selected AGN actually increases in galaxies in richer 
environments . This is comparable to results from the field by 
iLehmer et alJ ((2007), who studied X-ray emission from early- 



type galaxies in the Extended Chandra Deep Field South and 
found an average AGN fraction with hard X-ray luminosity 
above 10 41 erg s" 1 consistent with our measurement in clus- 
ters. For the broad X-ray band they find: /a (Mr < -20; Lx > 
10 41 ) - 7% and f A (M R < -20;L X > 10 42 ) = 2% (B. Lehmer 
2006, private communication). These fractions are remark- 
ably similar to the mean AGN fractions we measure in clus- 
ters and may suggest the X-ray luminous AGN fraction is 
not a strong func tion of environment, similar to the result of 
IBest et ail d2005l) for radio luminous AGN, although this does 
not expl ain the absence of X-ray luminous AGN in the group 
study of lShen et al.1 (120071) . 

AGN in clusters of galaxies have also garnered significant 
recent interest as an explanation of the absence of substan- 
tial cold gas at the centers of many clusters. Simple radiative 
cooling models predict that the intracluster medium (ICM) in 
the cores of many clusters should cool in less than a Hub- 
ble time, yet the predicted substantial reservoirs of cold gas 
are not observed. The presence of powerful radio gala xies at 
the ce nters of most predicted 'cool-core' clusters (e.g.. lBurnsl 
1990) suggests that while the cooling gas may provide fuel for 
the AGN, the AGN may also be reheating the cool gas. The 
substantial cavities in the hot ICM coincident with the lobes 
of these radio galaxies demonstrates t hat the AGN inject copi- 
ous amounts of en e rgy into the IC M (McNamara et aUbOOCt 
iFabian et alll2000h . iBirzan et all (12004 have shown that the 
energy necessary to reheat the ICM and prevent cooling is 
approximately consistent with the amount of energy neces- 
sary to create the large cavities observed in the ICM, although 
the mechanism by which the highly collimated AGN jets uni- 
formly heat the core ICM is still under active investigation. 
One possible resolution to the problem of unifor m heating is 
multip le AGN in the cluster core (~ 100 kpc). iNusser et alJ 
(2006) recently showed that gas cooling at the centers of clus- 
ters could condense onto any galaxies within the core, fuel ac- 
cretion onto the supermassive black holes, and produce more 
distributed heating of the ICM. A simple test of this scenario 
is to search for multiple or off-center AGN in the cores of 
clusters. It is also important to quantify the population and 
evolution of radio-bright AGN in the cores of clusters for fu- 
ture experiments that employ the Sunyaev-Zel'dovich effect 
to identify clusters be cause these source s could be a signifi- 
cant contaminant (e.g., Coble et al. 2003). 

In the present paper we employ the other cluster members 
identified in our multiwavelength survey to derive the spa- 
tial and kinematic distribution of the cluster AGN relative to 
the cluster galaxy population. We also use these data to de- 
rive velocity dispersions and study the amount of star forma- 
tion present in these clusters to determine if any global cluster 
properties correlate with the cluster AGN fraction. In the next 
section we provide a brief summary of the observations de- 
scribed in Paper I, followed by a derivation of the velocity 
dispersion, membership, and visible-wavelength properties of 
the cluster galaxies in § [3] The distribution of the AGN rela- 
tive to other cluster members is described in §|4]and the rela- 
tion between AGN fraction and the properties of the cluster is 
discussed in § [5] We present our conclusions in the last sec- 
tion. Throughout this paper we assume that the cosmological 
parameters are: {Hm, = (0.3, 0.7, 0.7) where Ho = 100/i 

km s _1 Mpc _1 . 

2. OBSERVATIONS 

The goal of this survey was to identify AGN in clusters 
of galaxies. This was accomplished with X-ray and visible- 
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TABLE 1 
Cluster Properties 



Cluster 


etc 


5c 


z 


Z1.Z2 


a 


'2(10 


Ps 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



Abell 3125 03:27:17.9 -53:29:37 0.0616 

Abell3128 03:30:43.8 -52:31:30 0.0595 

Abell644 08:17:25.6 -07:30:45 0.0701 

Abell2104 15:40:07.9 -03:18:16 0.1544 

Abelll689 13:11:29.5 -01:20:28 0.1867 

Abell2163 16:15:49.0 -06:08:41 0.2007 

MS1008 10:10:32.4 -12:39:53 0.3068 

AO 14 22:58:48.4 -34:48:08 0.3148 



0.0530,0.0700 475(94) 1.14 < 0.001 

0.0435,0.0755 906(74) 2.18 < 0.001 

0.0531,0.0871 952(382) 2.28 0.236 

0.1304,0.1783 1242(194) 2.85 0.203 

0.1392,0.2343 2402(357) 5.41 0.994 

0.1731,0.2236 1381 (324) 3.09 0.283 

0.2921,0.3215 1127 (153) 2.38 0.722 

0.2884,0.3412 2025 (217) 4.25 0.156 



NOTE. — Cluster sample and properties derived from the present study. Columns are: (1) Cluster name; (2 and 3) RA and DEC of the cluster center for epoch J2000; (4) redshift; 
(5) redshift range of cluster members; (6) velocity dispersion and uncertainty; (7) f200 in Mpc; (8) probability that the observed substructure is consistent with a random redistribution 
of the cluster radial velocities. The derivation of these quantities is discussed in Sections[3landl4l 
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FIG. 1. — R-band images of Abell 2104 (z = 0.154) and Abell 1689 (z = 0.187) with Chandra X-ray (black contours) and FIRST radio (white contours) data. 
Each image is 7' X 7' on a side, North is up and East is to the left. These images only show the central region of our data where the X-ray emission in prominent. 



wavelength images and spectroscopy of the eight clusters of 
galaxies listed in Table [TJ As AGN are quite luminous at X- 
ray wavelengths, AGN candidates were identified as galaxies 
with X-ray counterparts. These galaxies were then observed 
spectroscopically to determine if they were members of a 
given cluster. This X-ray selection technique was employed 
for two reasons. First, AGN are relatively rare in clusters 
of galaxies and it is more efficient to obtain complete spec- 
troscopic observations of all AGN candidates than to spec- 
troscopically observe all bright galaxies because the surface 
density of candidate cluster members with X-ray counterparts 
is substantially lower than the surface density of all candi- 
date cluster members. Secondly, X-ray selection is a rela- 
tively unbiased method of identifying AGN over a wide range 
of redshifts. Most AGN are low-luminosity (e.g., Lx < 10 43 
erg s" 1 ) and these AGN are difficult to identify spectroscopi- 
cally at even low-redshift (e.g., z ~ 0. 1) due to dilution by host 
galaxy starlight and/or obscuration. X-ray emission from an 
AGN has substantially higher contrast over other sources of 
X-ray emission to much lower accretion power, particularly 
when the spectroscopic aperture includes a substantial frac- 
tion of the total emission of the galaxy. X-ray emission is also 
much less affected by obscuration than visible-wavelength 



emission. 

The observing strategy was therefore to compare visible- 
wavelength and X-ray images of each cluster to identify po- 
tential cluster AGN for follow-up spectroscopy. As dis- 
cussed in greater detail in Paper I, the X-ray observations are 
archival Chandra observations of low-redshift clusters with 
sufficient sensitivity to identify low-luminosity AGN. We ob- 
tained visible-wavelength observations of the clusters with the 
2.5m du Pont telescope at Las Campanas Observatory with ei- 
ther the WFCCD or Tek5 CCD cameras. /?-band images of 
the centers of six of these clusters are shown in Figures [T] and 
|2] along with X-ray and radio contours. Comparable fig ures 
for A3 125 and A3 128 are presented in lRose et all (I2002I) . All 
X-ray counterparts brighter than R ~ 23 mag were then se- 
lected for multislit spectroscopic observations. As the sur- 
face density of X-ray counterparts was substantially less than 
the potential packing density of spectroscopic slits on a given 
mask, additional bright galaxies were observed to determine 
the spectroscopic properties of cluster members without X- 
ray emission, measure the fraction of cluster members with 
X-ray emission, and in several cases derive the velocity dis- 
persion of the cluster for the first time. 

Table [2] presents a catalog of all of our successful spectro- 
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scopic observations. As noted above, these targets were se- 
lected through a variety of algorithms. First, the highest pri- 
ority was assigned to all X-ray sources brighter than R ~ 23 
mag. Additional candidate cluster members were then tar- 
geted for spectroscopy based on brightness, color, and avail- 
ability of a spectroscopic slit. We successfully measured spec- 
troscopic redshifts for all sources brighter than a rest-frame 
absolute magnitude of Mr = -20 mag at each cluster redshift. 
For the nearby clusters, the apparent magnitude limits extend 
substantially fainter. In general, we are complete to R = 21 
mag for all of these fields and to R = 22 mag for all sources 
selected as X-ray counterparts. In addition, there are many 
fainter sources with emission-line redshifts. The exception is 
the Abell 2163 field, where the 90% completeness for spec- 
troscopy is approximately R = 20.5 mag for both selection cri- 
teria. With the exception of Abell 644 and Abell 2163, there 
are literature redshift data for all of these fields. We observed 
a number of these sources in each field and used these data to 
determine the typical redshift uncertainty of our observations 
is less than a z = 0.0005. 

3. CLUSTER PROPERTIES 

We use the membership data and the /?-band images and 
X-ray contours shown in Figures Q~]and[2]to either derive new 
centers for each cluster or confirm a previously-reported cen- 
ter from the literature. These cluster centers are listed in Ta- 
ble [U For most of the clusters, either a bright cluster mem- 
ber or the brightest cluster galaxy (BCG) is coincident with 
the approximate peak of the ICM. In these cases we adopt 
the coordinates of this galaxy for the cluster center, although 
for Abell 2163 the BCG candidate was not observed spectro- 
scopically. Only Abell 3 125 and Abell 3 128 do not have obvi- 
ous, bright cluster galaxies at their centers. Abell 3125 in fact 
does not have a diffuse ICM, while the ICM for Abell 3128 
is double-peaked. This pair of mergin g clusters have be en ex- 
tensively discussed and modeled by iRose et all (120021) . and 
we adopt their position for the center of Abell 3128, which 
is approximately midway between the peaks in the extended 
X-ray emission. For Abell 3125 we simply adopt the mean 
position of all of the confirmed membe rs, which is wit hin an 
arcminute of the position reported bv lAbell et al.l dl989). Fig- 
ure|3]shows the distribution of cluster members relative to the 
adopted centers along with the area subtended by the Chandra 
field of view. 

3.1. Velocity Dispersions and Membership 

We have used the spectroscopic data presented in Table [2] 
and literature membership data to determine the redshift, ve- 
locity dispersion, and redshift limits of each of these clusters. 
These values are presented in TableQ] We calculate these val- 
ues with an implementation of the center and scale estimators 
described in lBeers et al.l (|T990j). F° r each cluster we combine 
our redshift measurements and any from the literature and cal- 
culate an approximate redshift and velocity dispersion using 
the biweight estimator for center and scale. We then removed 
galaxies more than 5a from the center and iterated to deter- 
mine the mean redshift and velocity dispersion of each cluster. 
The adopted redshift range is ±5er of the cluster mean. The 
uncertainty in the velocity dispersion is the 90% confidence 
limit calculated with the jackknife of the biweight estimator. 
Figure |4] shows the galaxy velocity distribution for each clus- 
ter over the range ±3cr. 

As noted above, there are literature measurements of multi- 
ple cluster members for all of these clusters except Abell 644 



and Abell 2163. Ro se et ai1d200 2) performed a detailed study 
of the merging clusters Abell 3125 and Abell 3128 with spec- 
troscopic observations over a 2 degree field and mapped out 
the complex velocity structure of this cluster pair. Our ve- 
locity measurements are broadly consistent with their results, 
although given the kinematic complexity of this system the 
derived velocity dispersions are unlikely to be a good mea- 
sure of the cluster mass. As our obs ervations have a su bstan- 
tially smaller field of view than the IRose et aT] (|2002) study 
and map only a small fraction of the physical extent of these 
clusters, the present observations are not as well suited to 
map the velocity structure in detail. Our measurements for 
Abell 2104, Abell 1689, MS1008, and AC 114 are in good 
agreement with previous values. This is least surprising for 
Abell 2104, since we only have new redshif t s for six X-ray 
counterparts to add to the study of Liang et al. (2000). AC 1 14 
has the largest difference from the literature value. We mea- 
su re a 15% larger velocity di spersion than the earlier estimate 
by ICouch & Sharpies] dl987l) . although a substantial number 
of new members have been identified since that work. Our 
va lue for Abell 1 689 is comparable to the early measurement 
of iTeague et al.l d 1990b and larger than the recent value of 
Czoskg (120041) . who measure a w 2100 km s -1 based on a 
larger sample of over 500 (unpublished) members; however, 
the difference between these values has no impact on our re- 
sults below. 

3.2. Galaxy Colors and Spectral Properties 

In addition to changes in the mix of galaxy populations 
within clusters , such as the morphology-density relation 
dDresslerlll98 0'). the mix of galaxy populations can also vary 
between cluster s. One clear example of this is the Butcher- 
Oemler effect ( Butcher & Oemlerj 1 1978b . the tendency for 
clusters to have lar ger blue galaxy fr a ctions at higher red- 
shift. Specifically, iButcher & Oemlerj d 19781) quantified the 
blue galaxy fraction in a cluster as the fraction of galax- 
ies more than 0.2 mag bluer in rest-frame B - V color. The 
color-magnitude diagrams presented in Figure [5] indicate that 
the confirmed cluster members (triangles) have a substantial 
variation in observed B-R color. We calculate the Butcher- 
Oemler (BO) galaxy fraction for each cluster as the fraction 
of all known members more luminous than Mr < -20 mag in 
our photometric data set. The number of BO galaxies and BO 
fraction are listed in Table [3] The highest fractions are mea- 
sured for clusters with substantial starforming galaxy popula- 
tions (e.g., A1689, AO 14). While the BO fraction is likely 
sensitive to potential biases such as the fraction of the cluster 
surveyed (e.g., fraction of the virial radius) and the limiting 
magnitude of the spectroscopy, the relative fractions provide 
a measure of the relative star formation activity from cluster 
to cluster in this sample. The difference between our mea- 
su red BO fraction of 1 8% fo r AC1 14 and the 29% measured 
bv lCouch & Sharpies! (1 19871) provides a reasonable gauge of 
how sensitive this quantity is to measurement details. 

We have also used our spectroscopic data for all of the clus- 
ters except Abell 2104 to determine the emission-line galaxy 
fraction for galaxies more luminous than Mr < -20 mag. This 
provides a second measure of the frequency of star formation 
in cluster members, although this quantity is less reliable than 
the BO fraction because emission-lines may be more diffi- 
cult to detect in cluster members at higher-redshift because 
the typical signal-to-noise ratio of the spectroscopy is lower. 
We simply classify any galaxy with an emission line as an 
emission-line galaxy, regardless of the strength of the line. In 
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TABLE 2 
Spectroscopic Catalog 



RA 


DEC 


z 


Template 


Mask 


Select 


R 


B-R 


V-R R-I 


Lit ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) 


(11) 


J ,JLi 1 .'+~> . J *J 


-51-74-D7 6 


0599 




a3125a 


p 


1 6 RQ (0 Oil 


1 65 (0 fi51 


5Q (0 051 




1-97-50 2? 


-51-24-54 7 


0.5402 


f£ 


a3125a 


p 


70 67 (0 Oil 


1 72 (0 1 01 


94 (0 061 






-51-21-52 1 


0.0626 


A 


a3125a 


p 


1 7 17 10 Oil 


1 66 (0 051 


59 (0 051 




3:27:22.03 


-53:25:57.6 


0.3950 


E 


a3125a 


P 


21.14(0.04) 




0.87 (0.11) 




3:27:52.60 


-53:24:08.4 


0.0614 


A 


a3125a 


p 


15.93 (0.03) 


1.77 (0.05) 


0.66 (0.05) 


2MASX J03275262-5324079 


3:27:46.86 


-53:22:58.7 


0.2581 


E 


a3125a 


P 


20.16(0.03) 


1.62 (0.06) 


0.64 (0.05) 




3:27:46.88 


-53:22:08.8 


0.3937 


E 


a3125a 


P 


20.29 (0.03) 


1.66 (0.06) 


0.69 (0.05) 




3:27:56.19 


-53:23:18.4 


0.4374 


E 


a3125a 


P 


20.88 (0.03) 


1.40 (0.08) 


0.68 (0.05) 




3:27:38.88 


-53:26:04.5 


0.597 


Q 


a3125a 


X 


19.91 (0.03) 


0.35 (0.05) 


0.14(0.05) 




3:27:54.32 


-53:21:51.1 


0.84 


Q 


a3125a 


X 


20.80 (0.03) 


0.64 (0.05) 


0.28 (0.05) 




3:26:55.66 


-53:31:58.8 


0.3192 


A 


a3125b 


P 


19.18 (0.03) 


2.71 (0.07) 


0.97 (0.05) 




3:26:54.48 


-53:31:47.9 


0.3101 


E 


a3125b 


P 










3:27:00.96 


-53:31:36.6 


0.0630 


E 


a3125b 


P 


19.91 (0.03) 


1.09 (0.05) 


0.38 (0.05) 




3:27:30.01 


-53:31:04.4 


0.2613 


A 


a3125b 


P 


19.72 (0.03) 


1.55 (0.05) 


0.56 (0.05) 




3:27:22.70 


-53:30:32.9 


0.5374 


A 


a3125b 


P 


20.75 (0.03) 


1.29 (0.07) 


0.84 (0.06) 





NOTE. — Spectroscopic catalog of all sources with successful redshift measurements. Columns are: (1 and 2) object RA and DEC in J2000; (3) redshift; (4) best template 
match as emission-line galaxy (E), absorption-line galaxy (A), quasar (Q), or Galactic star (S); (5) mask identification; (6) target selection as either an X-ray source or from 
R-band photometry; (7) R-band magnitude and uncertainty; (8 to 10) B-R, V-R, and R-I color and uncertainty; (11) literature identification. The complete version of this 
table is in the electronic edition of the Journal. The printed edition contains only a sample. 



TABLE 3 

Cluster Galaxy Population Fractions 



Cluster 


N B o 


/bo 


N e 


fe 


(1) 


(2) 


(3) 


(4) 


(5) 


A3125 


2 


0.11 


4 


0.44 


A3128 


5 


0.07 


7 


0.21 


A644 


2 


0.13 


1 


0.06 


A2104 


3 


0.09 






A 1689 


31 


0.25 


8 


0.24 


A2163 


3 


0.12 


3 


0.11 


MS 1008 


18 


0.31 


4 


0.16 


AO 14 


25 


0.26 


2 


0.18 



NOTE. — The fraction of Butcher-Oemler and emission-line galaxies in 
each cluster that are brighter than a fixed rest-frame absolute magnitude of 
Mr = —20 mag. For each cluster listed in column 1 we list the total num- 
ber of Butcher-Oemler galaxies that are known cluster members in column 2 
and their fraction of the total cluster galaxy population above that luminos- 
ity in column 3. We present the total number of emission-line galaxies and 
the emission-line galaxy fraction in columns 4 and 5, although this quantity 
is only calculated from galaxies we observed spectroscopically. We do not 
have sufficient data to calculate this quantity for A2104. See § !3.2l for further 
details. 

the vast majority of cases, the observed emission-line is the 
[Oil] A3727 doublet and nearly all of the emission-line galax- 
ies appear to be starforming galaxies. The few exceptions 
are the small number of X-ray-identified AGN with visible- 
wavelength emission-line signatures. We do not spectroscop- 
ically identify any AGN that do not have an X-ray counter- 
part. Both the BO fraction and emission-line galaxy fraction 
for these clusters are compared with the AGN fraction in Sec- 
tion|5]below. 

4. DISTRIBUTION OF CLUSTER AGN 

Major mergers between gas-rich galaxies remain the stan- 
dard paradigm for triggering high-luminosity quasars, al- 
though there is little direct evidence that such spectacular 
events are resp onsible for triggering lower-lum inosity AGN 
as well ("e.g.. Ide Robertis et alJfl998t ISchmitll2001h . How- 
ever, in the standard picture of merger-driven fueling, the 
AGN luminosity dim inishes gradually over up to a Gyr (e.g., 
iHopkins et alj |2005). Little evidence of a past merger may 





TABLE 4 






KS Test Results 




Distribution Type 


Sample Limits 


KS 


(1) 


(2) 


(3) 



Velocity distribution 


Lx 


> 


10 41 




0.683 




Lx 


> 


10 42 




0.955 




Lx 


> 


10 41 - 


-A3125/8 


0.388 




Lx 


> 


10 42 - 


-A3 125/8 


0.967 


Radial [Mpc] 


Lx 


> 


10 41 




0.743 




Lx 


> 


10 42 




0.069 




Lx 


> 


10 41 - 


-A3125/8 


0.950 




Lx 


> 


10 42 - 


-A3125/8 


0.128 


Radial [r 20 o] 


Lx 


> 


10 41 




0.989 




Lx 


> 


10 42 




0.028 




Lx 


> 


10 41 - 


-A3125/8 


0.513 




Lx 


> 


10 42 - 


-A3 125/8 


0.051 


Substructure 


Lx 


> 


10 41 




0.616 




Lx 


> 


10 42 




0.606 




Lx 


> 


10 41 - 


-A3125/8 


0.302 




Lx 


> 


10 42 - 


-A3 125/8 


0.649 



NOTE. — Results of the Kolmogorov-Smirnoff tests. Columns (11 and (2) 
describe the input sample of X-ray sources; in all cases the comparison sam- 
ple comprises the cluster members without X-ray emission. Column (3) lists 
the probability that the X-ray sources are drawn from a different distribution. 
Sample limits marked — A3125/8 did not include these two clusters. See §|4*1 
for further details. 



therefore be present during a more extended, low-luminosity 
phase. 

As outlined above, major mergers between gas-rich galax- 
ies should be substantially rarer in clusters of galaxies than in 
the field because the velocity dispersion is too high for major 
mergers and few galaxies are gas rich. The most favorable 
region for major mergers to still occur in clusters is in the 
outskirts, which contain many galaxies falling into the clus- 
ter for the first time. These galaxies are primarily spirals and 
therefore are richer in cold gas. They are also often bound in 
low-velocity groups with velocity dispersions of a few hun- 
dred kilometers per second and consequently have a higher 
merger probability. 

Previous studies have shown that there is more activity in 
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FIG. 2.— R-band images of Abell 644 (upper left, z = 0.070), Abell 2163 (upper right, z = 0.201), MS1008 (lower left, z = 0.307), and AC114 (/ovver rig/if, 
z = 0.315) along with Chandra X-ray (black contours) and NVSS radio (white contours) data. Each image is 7' X 7' on a side, North is up and East is to the left. 
As in Figure[T] these images only show the central region of the cluster where the X-ray emission is prominent. 



the outskirts of clusters, such as a la rger fraction of emission- 
line and poststa rburst galaxies (e.gJCouch & Sharplesj[l98Zt 
Fabri cant etaOl991t iFisher et al.ll 19981) . although these stud- 
ies did not specifically address the AGN distr i bution . In a 
study of 10 low-redshift clusters, Dress leFet al.l (fl999) found 
that emission-line galaxies had a higher velocity dispersion 
than the more passive galaxies. The higher line-of-sight ve- 
locity dispersion implies that the emission-line galaxies are 
not yet virialized and remain on primarily radial orbits in the 
cluster potential. 

Clusters of galaxies offer a unique opportunity to test if low- 
luminosity AGN are primarily the result of the major mergers 
of gas-rich galaxies because these mergers are most likely to 
occur in the infalling population and the cluster crossing time 
is comparab le to the predi cted lifetime of the low-luminosity 
phase (e.g.. iMartinil 120041) . Low-luminosity AGN triggered 
in infalling galaxies will therefore remain on primarily radial 
orbits after a Gyr, long after evidence of their violent past has 
faded. In the next subsections we investigate the kinematic 
and radial distribution of the cluster AGN to determine if they 



may have recently entered the cluster potential. 

4.1. Velocity Distribution 

In Figure |4] we show the galaxy velocity distribution for 
each cluster and mark the locations of all of the X-ray sources 
with arrows and a dashed histogram. While there are a num- 
ber of X-ray sources several standard deviations from the 
mean cluster redshift, there are an insufficient number of 
sources per cluster to state if these outliers represent a sig- 
nificant fraction of the population. 

To increase the AGN sample size for statistical analysis, 
we have combined the velocity distributions for all eight clus- 
ters. Because there is a substantial variation in the number of 
galaxies, number of X-ray sources, and velocity dispersions, 
we combined these clusters by first calculating the absolute 
value of the velocity offset of each galaxy from the cluster 
mean with respect to the cluster velocity dispersion. All eight 
clusters were then summed in these normalized coordinates 
and in Figure [6] we show the cumulative distribution of the 
X-ray sources relative to the cluster galaxies without X-ray 
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FIG. 3. — Positions of cluster members relative to the centers listed in Tablefl] Members are marked with different symbols if their positions are from the 
literature (small, black circles), identified in the present study (large, red circles), and X-ray sources (blue crosses). Also shown are the field of view of the 
Chandra observations (solid, black boxes) and circles of radius 0.5 Mpc (smaller, dotted circle), 1 Mpc (larger, dotted circle), and r2oo (large, dashed circle), if 
the later radius falls within the field of view shown. The axes on all panels are offsets in arcminutes from the cluster center. 



emission. We plot separate velocity distributions for all clus- 
ter AGN more luminous than Lx > 10 41 erg s" 1 and all more 
luminous thanLx > 10 42 erg s" 1 . 

This figure demonstrates that the velocity distribution of 
cluster AGN is essentially identical to the other cluster mem- 
bers. A Kolmogorov-Smirnov (KS) test confirms this to be 
the case, specifically that there is a 68% and 96% probability 
that the L x > 10 41 erg s _1 and L x > 10 42 erg s" 1 AGN, re- 
spectively, are drawn from the same parent population as the 
cluster members without X-ray counterparts (see Table|4|i. We 
therefore do not find evidence that the cluster AGN are pref- 
erentially distributed on radial orbits, even if we exclude the 
merging clusters A3 125/8, unlike the ca se for emission-line 
galaxies in clusters (Dres sier et al.ll 1999b . This may indicate 
that nuclear activity in cluster galaxies can remain, or be reac- 
tivated, after the parent population has virialized in the cluster 
potential. 

To test the sensitivity of our data to the known differences 
between cluster galaxy populations, we computed similar dis- 



tributions between Butcher-Oemler and non-Butcher-Oemler 
galaxies as well as between emission-line and absorption- 
line galaxies. This analysis showed that the Butcher-Oemler 
galaxies have systematically larger velocity dispersions than 
the redder cluster members, as expected for a population on 
more radial orbits. The difference between emission-line and 
absorption-line cluster members was not statistically signifi- 
cant, although as noted previously the emission-line classifi- 
cation is more susceptible to variations in the average signal- 
to-noise ratio than the color classification. The emission and 
absorption-line samples used here are also only based on our 
spectroscopy, and not all known members, and therefore the 
sample is smaller than that used for the Butcher-Oemler and 
red galaxy comparison. 

4.2. Radial Distribution 

To investigate the radial distribution of the AGN relative 
to the cluster center, we have similarly combined all of the 
cluster sources shown in Figure [3] and plot the cumulative 
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FIG. 4. — Velocity distributions for each cluster. Each panel displays the distribution of all galaxies within 3cr of the cluster mean redshift (black, solid line), 
the subset of new members (red, dotted line), and the X-ray sources (blue, dashed line) for a cluster field. The redshifts of individual X-ray sources are also 
marked (blue arrows). The bottom abscissa plots the observed redshift, while the top abscissa plots the rest-frame velocity offset from the cluster mean. 

radial distribution of the AGN and inactive cluster members 
in Figure [7] We have calculated this distribution as a func- 
tion of projected physical distance from the cluster center in 
Mpc and as a function of projected fraction of r2oo, where 
r2oo is the physical radius within which the mean density of a 
virialized cluster of the measured velocity dispersion exceeds 
the critical density at that redshift b y a factor of 200 (e.g., 
ICarlberg et aHll997tlTreu et al.ll2003l) . This figure shows that 
the Lx > 10 41 erg s" 1 cluster AGN and the other cluster mem- 
bers have a similar radial distribution, but above Lx > 10 42 
erg s the cluster AGN are more strongly concentrated. A 
KS test confirms the visual impression of both panels of Fig- 
ure|7] There are formally only 7% and 3% probabilities that 
the Lx > 10 42 erg s" 1 AGN are drawn from the same distribu- 
tion in metric and virial radius, respectively (see Table|4]i. The 
probabilities are less than a factor of 2 higher when we cluster 
the merging clusters A3 125/8. This result confirms the previ- 
ous, purely statistical e vidence for a concentration of AGN in 
the ce nters of clusters dRuderm an & Ebeling 20QH iDowsettl 
120051) with a spectroscopically-confirmed sample. In a recent 



paper [On & Mohrl (|2007) found evidence that radio sources 
in clusters are also more centrally concentrated than typical 
cluster galaxies, however we intriguingly do not find any ev- 
idence that the X-ray and radio AGN are in the same host 
galaxies. Compariso n of our X-ray selecte d AGN and those 
in the radio study of Mor rison et al.l (120031) indicates that the 
two samples are nearly disjoint (see also the radio and X-ray 
contours in Figure [TJ. 

The summed radial distribution of cluster members is more 
sensitive to several potential biases than the summed veloc- 
ity distribution. These biases include the physical area of the 
cluster surveyed for X-ray and other members and the relative 
physical sizes of the clusters. The first is important because 
the distributions will artificially appear different if the field of 
the X-ray observations is larger or smaller than the field of 
view surveyed to obtain membership information for inactive 
galaxies. We have addressed this source of systematic error 
by only including our spectroscopically confirmed members 
in Figure [7] rather than all available members. Our spectro- 
scopic sample is a suitable comparison sample because the 
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R magnitude 

FIG. 5. — Color-magnitude diagrams for the eight cluster fields. All sources in our R— band catalog with measured B-R colors and R-band stellarity less than 
0.9 are shown. Confirmed members from Table[2]trace the color-magnitude relation of the passively-evolved galaxies in each cluster (red, open triangles). Other 
members from the literature (small triangles) and X-ray sources (blue crosses) are also marked. The vertical, dashed line marks the apparent magnitude that 
corresponds to a rest-frame Mr = -20 mag for each cluster. 



X-ray and other cluster galaxy candidates were observed with 
the same multislit masks and the AGN and inactive galaxies 
in each cluster were targeted over the same radial distribution 
from the cluster center. We further refined the input catalog 
by only including confirmed cluster members that fall within 
the field of view of the Chandra observations. 

The second potential bias is the different sizes of the clus- 
ters, both their physical size and projected angular extent rel- 
ative to the Chandra field of view. For example, the ACIS-I 
camera on Chandra only encompasses the central ~ 0.5 Mpc 
of Abell 3125 and Abell 3128, while it encompasses up to 
~ 3 Mpc from the center of MS 1008 (see Figure[3]l. The red- 
shift range of the sample therefore causes the higher-redshift 
clusters to dominate the distribution at the largest distances 
from the cluster center. Similarly, only the low-mass cluster 
A3 125 and the higher-redshift cluster MS 1008 provide sub- 
stantial information about the distribution of cluster galaxies 
at more than half of r2oo- The virial scaling does not change 
the fact that the Chandra field of view samples a a different 



fraction of the physical extent of different clusters. 

A final potential bias on the radial distribution of X-ray 
sources is the detection efficiency of X-ray sources as a func- 
tion of distance from the cluster center. The detection ef- 
ficiency may be a function of clustercentric distance due to 
two effects: the higher background near the cluster center due 
to the ICM (see Figures Q] and [2]) and the degradation of the 
Chandra point-spread function (PSF) at larger off-axis angles 
may make it more difficult to detect sources. While the effect 
of the bright ICM would decrease the probability of detecting 
an X-ray source near the cluster core and therefore is opposite 
our observed trend, the larger Chandra PSF at larger off-axis 
angles could mimic the observed trend for Lx > 10 42 erg s" 1 
AGN. 

To quantify both of these potential biases on the apparent 
radial distribution, we performed a series of simulations with 
the MARX 1 package to quantify our detection efficiency as 

1 http://space.mit.edu/CXC/MARX 
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FIG. 6. — Cumulative, normalized velocity distribution of cluster X-ray 
sources relative to each cluster's mean redshift. The distribution of the AGN 
more luminous than Lx > 10 4 ' erg s~ 1 (blue, dotted line), more luminous than 
Lx > 10 42 erg s (blue, dashed line), and the inactive galaxies (red, solid line) 
are statistically indistinguishable. The absolute value of the velocity offset 
of each galaxy from the clus ter mean has been normalized to the cluster's 
velocity dispersion. See 34.1l for further details. 

a function of source luminosity and clustercentric distance. 
MARX is well-suited to this task because it can generate a 
PSF that accounts for the off-axis angle, aspect solution, and 
spectral shape. The detection efficiency relative to the ob- 
served radial distribution of all members with Mr < -20 mag 
and no X-ray detection is the most relevant quantity because 
our comparison is directly with these galaxies. We therefore 
used MARX to add an artificial X-ray point source at the 
location of every cluster member without an X-ray source. 
This source was assigned a T = 1 .7 spectra shape and the total 
counts corresponding to a given luminosity at the redshift of 
the particular cluster. The fraction of these sourc es recovered 
with th e same WAVDETECT parameters used in Mart ini et alj 
(2006) then provides a measurement of the detection effi- 
ciency. We performed this exercise for point sources with 
total counts corresponding to a source with L x = 10 41 and 
Lx = 10 42 erg s" 1 for each cluster and then calculated the de- 
tection efficiency as a function of clustercentric distance in 
both Mpc and fraction of the virial radius. From this approach 
we find that the detection efficiency varies by less than 20% 
out to 1.5Mpc and 0.5r2oo, after which the results begin to be 
compromised by small number statistics (few cluster mem- 
bers). This is true for both 10 41 and 10 42 erg s" 1 sources and 
we therefore conclude that the observed trend with radius is 
real and not an artifact of either bias. We note further that 
this approach is somewhat pessimistic as it assumes that all of 
the point sources are at the luminosity limit, rather that at or 
above it, but it nevertheless provides a relative measure of the 
detection pr obab ility. 

As in § 14.11 we compared the radial distributions of 
photometrically (Butcher-Oemler or red) and spectroscopi- 
cally (emission-line or absorption-line) defined subsamples 
to test the sensitivity of these data to other known rela- 
tions between cluster populations. This analysis showed 
that the red galaxies are more centrally-concentrated than the 
Butcher-Oemler population, as expected. We also found no 
statistically-significant difference between the emission-line 
and absorption-line populations, which we again attribute to 



the smaller size of the subsamples input to this analysis. 
4.3. Association with Cluster Substructure 

If cluster AGN are preferentially associated with infalling 
groups of galaxies, or other structures within clusters, then 
they ma y trace local deviations fro m the mean cluster kine- 
matics. Dressier & Shec tmanl (119881) developed a substruc- 
ture test to identify deviations in the cluster velocity field, 
where the substructure parameter S is defined as 

S 2 = ( 1 1 /a 2 ) [(v local — v) 2 + (a local -a) 2 ] ( 1 ) 

and the local values of the radial velocity v\ oca i and velocity 
dispersion a\ oca i are calculated from the radial velocities of 
each galaxy and its 10 nearest neighbors. Larger values of S 
correspond to larger local deviations from the cluster kinemat- 
ics. Figure [8] plots all of the cluster members shown in Fig- 
ure |3] with open circles scaled by this s ubstructure measure. 
One m inor change from the method of Dre ssier & Shectmanl 
( 1988) is that we calculate these quantities with the biweight 
estimator described previously, rather than assume a Gaussian 
distribution. 

This analysis indicates that Abell 3125 and Abell 3128 
have th e most substantial substructure. Dres sier & Shectmanl 
( 1988) note that the sum of the 5 values for all cluster mem- 
bers is a measure of how much substructure is present and 
a simple way to characterize this value for a given cluster is 
via Monte Carlo simulations with random reshuffling of the 
measured velocities among the measured positions. We have 
generated 10,000 realizations with random redistributions of 
the measured velocities among the measured galaxy positions 
and calculated the probability of obtaining the measured sum 
of the 5 values or larger. These values are listed in Table[T]and 
indicate that Abell 3125 and Abell 3128 indeed have signifi- 
cant substructure, which is not surprising for a merging pair of 
clusters, while the remaining clusters do not have substantial 
substructure. 

The substructure parameter can also be used to character- 
ize the extent that a particular cluster population is associated 
with substructure. We have calculated the distribution of S 
values for the X-ray sources in all eight clusters relative to 
the other cluster members and plot this distribution in Fig- 
ure [9] There is good agreement between the distribution of 
the Lx > 10 41 erg s" 1 , Lx > 10 42 erg s" 1 , and cluster members 
without X-ray counterparts, and this indicates that the cluster 
AGN are not more likely to be associated with substructure 
in clusters, at least on the scales this parameter is sensitive to 
(a galaxy and its ten nearest neighbors). As Abell 3125 and 
Abell 3128 have substantial substructure, we repeated this ex- 
ercise without these two cluster and still obtained consistent 
distributions. Results from KS tests of all of these distribu- 
tions are listed in Table |4] 

5. RELATION OF AGN TO CLUSTER PROPERTIES 

In the previous two sections we have analyzed the proper- 
ties of the clusters and the distribution of the AGN within the 
clusters. In the present section we discuss the completeness 
of our survey of each cluster and investigate potential corre- 
lations between the AGN fraction and the properties of their 
host cluster. 

5.1. Completeness 

In Paper I we calculated the AGN fraction in clusters of 
galaxies and defined the AGN fraction as the fraction of galax- 
ies more luminous than a rest-frame Mr = -20 mag with X- 
ray counterparts more luminous than Lx = 10 41 erg s" 1 in the 
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FIG. 7. — Cumulative, normalized radial distributions relative to each cluster's center in projected Mpc (left) and in units of r2oo (right). The distribution of 
the AGN more luminous than Lx > 10 42 erg s (blue, dashed line) is substantially more centrally concentrated than the inactive galaxies (red, solid line) in both 
radial distributions. The distribution of the AGN more luminous than Lx > 10 4 ' erg & {blue, dotted line) is comparable to the inactive galaxies. The ra dial 
distribution is calculated from the projected distance of each galaxy from the cluster center. Larger separations are sampled by the more distant clusters. See 34.21 
for further details. 



TABLE 5 

Completeness and AGN Fractions 



Cluster 
Name 



All 



-Nx «• 
> 10 41 



-20 mag) - 



>10 42 



Membership 

Confirmed Estimated 



All 



> 10 41 



•/a- 



10 42 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Abell3125 


6 


5 


4 





18 


28 


0.179 


0.143 


0.000 


Abell3128 


10 


10 


7 


1 


54 


63 


0.159 


0.111 


0.016 


Abell644 


2 


2 


2 


1 


15 


75 


0.027 


0.027 


0.013 


Abell2104 


6 


4 


4 


1 


33 


62 


0.065 


0.065 


0.016 


Abelll689 


2 


2 


2 


2 


110 


337 


0.006 


0.006 


0.006 


Abell2163 


3 


2 


2 


1 


25 


262 


0.008 


0.008 


0.004 


MS 1008 


5 


4 


4 


3 


75 


346 


0.012 


0.012 


0.009 


AC114 


6 


5 


5 


3 


95 


204 


0.025 


0.025 


0.015 


Average 














0.060 


0.049 


0.010 


Sum 


40 


35 


30 


12 


425 


1377 


0.025 


0.022 


0.009 



NOTE. — Completeness estimates and AGN fractions for galaxies more luminous than Mr = -20 mag and different X-ray luminosity cuts. Columns are: (1) Cluster name; (2) 
Number of X-ray sources in the cluster regardless of the host galaxy luminosity; (3) Number of X-ray sources in hosts more luminous than Mr = -20 mag; (4) Number of X-ray 
sources more luminous than Lx = 10 41 erg s -1 in the broad X-ray band; (5) Number of X-ray sources more luminous than Lx = 10 42 erg s -1 ; (6) Number of spectroscopically 
confirmed cluster members within the Chandra field of view more luminous than Mr = -20 mag; (7) Estimate of the total number of members more luminous than Mr = —20 mag; 
(8) AGN fraction if all X-ray sources are AGN; (9) AGN fraction if all X-ray sources more luminous than Lx = 10 41 erg s -1 are AGN; (10) same as column (9) for X-ray sources 
more luminous than Lx = 10 42 erg s -1 . The final two rows present the AGN fraction calculated from the average of the eight clusters and the average of the sum of the eight clusters. 
See § |5.1| for further details. 



broad X-ray band, or f A (M R < -2Q;L X > 10 41 ) = 5 ± 1.5%. 
The quoted uncertainty in this value was only a Poisson esti- 
mate and did not take into account potential systematic un- 
certainties, such as the estimate of the total cluster galaxy 
population. The determination of the numerator in the AGN 
fraction requires redshift measurements for all X-ray counter- 
parts brighter than the apparent magnitude of a cluster mem- 
ber with this absolute magnitude. As we discussed in §|2] our 
spectra are complete to the appropriate apparent magnitude 
limit for each cluster. Measurement of the denominator re- 
quires a much larger number of spectra because most galaxies 
to the requisite apparent magnitude limit in the Chandra field 
of view are not detected X-ray sources. For all but the lowest- 



redshift clusters, these data are substantially incomplete and 
we have to sum the known members with an estimate of the 
additional cluster members without redshifts to estimate the 
total number of cluster members brighter than M« = -20 mag. 
As noted in Paper I, we can quantify this completeness with 
photometric observations of known cluster members, known 
nonmembers, and galaxies without spectra. Here we estimate 
the total number of cluster members without spectra by cal- 
culating the fraction of galaxies with spectra that are clus- 
ter members as a function of B-R color, rather than adopt a 
color-independent completeness estimate as in Paper I. This 
color-dependent completeness estimate takes into account the 
observation that galaxies with similar colors to observed clus- 
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ter members are more likely to be cluster members than those 
with dramatically different colors, but still accounts for a po- 
tential population of bluer or Butcher-Oemler galaxies if such 
a population has been spectroscopically confirmed. A disad- 
vantage of this approach is that the completeness estimate is 
sensitive to the number of members and nonmembers with 
spectra at a given color. We have increased the number of 
galaxies with spectra at a given color by also including known 
members and nonmembers from the literature for which we 
have photometry, although they fall outside of the Chandra 
field of view. However, these additional data are only used to 
provide an improved completeness fraction as a function of 
color; the total number of cluster members is estimated only 
from galaxies that fall within the Chandra field of view. 

In Table [5] we use the results of this completeness calcula- 
tion to obtain an estimate of the total number of cluster mem- 
bers more luminous than Mr = -20 mag. For each cluster 
we list both the number of known members from our data 
and the literature that fall within the Chandra field of view 
and our estimate of the total cluster galaxy population. We 



then calculate the AGN fraction for these clusters with three 
different cuts in X-ray luminosity: (1) the fraction of AGN 
if all of the X-ray sources in galaxies with Mr < -20 mag 
identified in Paper I are AGN; (2) the AGN fraction for only 
those X-ray sources above a luminosity of Lx = 10 41 erg s , 
and (3) the AGN fraction for those above Lx = 10 42 erg s" 1 . 
The second of these three limits, namely Mr < -20 mag and 
Lx > 10 41 erg s" 1 corresponds to the AGN fraction quoted in 
Paper I. Here we derive f A (M R < -20;L X > 10 41 ) = 5%, in 
excellent agreement with our calculation from Paper I. This 
correction provides a better estimate of the membership frac- 
tion for galaxies observed to be redder and bluer than typical, 
known members. The X-ray luminosity threshold was chosen 
previously because it is rare for other potential sources of X- 
ray emission, most notably LMXBs and hot, gaseous halos, 
to produce such high X-ray luminosities. However, as noted 
in Paper I, the disadvantage of this threshold is that the X-ray 
observations of Abell 1689 and MS 1008 are not quite sensi- 
tive to this luminosity limit and therefore some AGN may be 
missed and the AGN fraction underestimated. We have also 
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FIG. 9. — Cumulative distribution of the substructure parameter S for X-ray 
sources more luminous than Lx > 10 41 erg s^'iblue, dotted line), more lumi- 
nous than Lx > 10 42 erg s~ 1 (blue, dashed line), and inactive cluster members 
(red, solid line). Larger S corresponds to a larger local deviation from the 
mean cluster kinematics. 

calculated the AGN fraction for a factor of ten higher X-ray 
luminosity and find f A (M R < -20;L X > 10 42 ) = 1%. All of 
our X-ray observations are sensitive to such luminous sources. 
Galaxies that can produce such luminous X-ray emission by 
mechanisms other than black hole accretion are also much 
rarer than AGN. 

The AGN fraction may also be a function of host galaxy lu- 
minosity, in add ition to X-ray lumin osity. In our original work 
on Abell 2104 dMartini et al. 2002), we noted that two of the 
three brightest cluster members were AGN. From inspection 
of Figure [5] a similar tendency is apparent at the bright end of 
the color-magnitude relation in Abell 3125, Abell 3128, and 
Abell 2163. In fact, approximately half (19) of the AGN more 
luminous than Lx > 10 41 erg s" 1 are in host galaxies more lu- 
minous than Mr < -21.3 mag and for these galaxies we find 
f A (M R < -2\3;L X > 10 41 ) = 9.8%. This absolute magnitude 
li mit was chosen t o correspond to the luminosity threshold 
of ISun et all (120071) . who find nine AGN with L x > 10 41 erg 
s" 1 in 163 galaxies in a sample of z = 0.01-0.05 clusters, or 
f A (M R < -2\3;L X > 10 41 ) = 5%. While this value is a factor 
of two below our value, the differences may be due to varia- 
tions in the AGN fraction between cluster samples (e.g. dif- 
ferent redshift distributions) and we discuss this point further 

in en 

An important caveat to these estimates is that they are es- 
timates of the AGN fraction within the field of view of the 
Chandra ACIS-I or ACIS-S detectors. These detectors have 
an approximate field of view of 16.9' x 16.9' and 8.3' x 8.3', 
respectively. These two angular sizes, combined with the red- 
shift distribution of our sample and the range in cluster sizes 
(e.g., see the r2oo value in Table [TJ, correspond to substan- 
tial variations in the fraction of each cluster galaxy popula- 
tion surveyed. This point is illustrated with Figure [3] which 
includes circles that mark the r2oo radius. The Chandra field 
of view approximately covers out to the r2oo radius for only 
MS 1008, while for many of the remaining clusters it is even 
larger than the 25' square panels shown for each cluster. Our 
measurements of the AGN fraction for each cluster are there- 
fore necessarily measurements dominated by the center of the 



cluster for nearly all cases. However, the absence of a strong 
radial trend in the AGN distribution for the more common 
L x > 10 41 erg s" 1 sources shown in Figure [7] suggests that the 
ratio of AGN to inactive galaxies may not be a strong function 
of position within the cluster. The AGN fraction may there- 
fore be relatively well determined, even if the census of the 
total cluster AGN population is incomplete. We also note that 
for several clusters our ground-based images do not encom- 
pass the entire Chandra field, and in fact some of the X-ray 
sources identified in A3 125 and A3 128 are from literature po- 
sitions and redshifts, but the measurement of the AGN frac- 
tion remains relatively robust in the absence of strong radial 
gradients. 

We have assumed that each cluster of galaxies has the same, 
intrinsic AGN fraction (and the same X-ray luminosity func- 
tion) and treated each cluster as an independent measure of 
the AGN fraction. This yields an average fraction of 5% for 
L x > 10 41 erg s" 1 and 1% for L x > 10 42 erg s _1 . If every 
cluster galaxy had an equal probability of hosting an AGN, 
an alternate way to calculate the AGN fraction in clusters of 
galaxies would be to treat each cluster galaxy as an indepen- 
dent measure of the AGN fraction, rather than each cluster. 
This yields an AGN fraction of 2.2% for L x > 10 41 erg s _1 
and 0.9% for L x > 10 42 erg s _1 . The value for L x > 10 41 
erg s" 1 sources is substantially lower than the average of the 
cluster AGN fractions and provides evidence that the AGN 
fraction may vary from cluster to cluster. 

5.2. Evidence for Intrinsic Variation 

We now investigate variations in the AGN fraction from 
cluster to cluster and consider if these variations are consistent 
with Poisson fluctuations and systematic errors in the clus- 
ter membership. If we ignore potential systematic errors in 
the completeness correction, the AGN fractions in these eight 
clusters are inconsistent with Poisson fluctuations with greater 
than 95% confidence for Lx > 10 41 erg s" 1 ; there is no evi- 
dence for variations in the AGN fraction for the much smaller 
population with Lx > 10 42 erg s" 1 . The dominant source of 
systematic error in the AGN fraction is the estimate for the 
number of cluster galaxies in the Chandra field of view. This 
correction is more important for clusters that have a larger 
number of galaxies within the field of view of the Chandra 
observation and those clusters that have less membership data. 
The ratio of the estimated to confirmed cluster galaxy popu- 
lations brighter than M R < -20 mag in Table [5] provides the 
size of the correction and varies from ~ 1 for Abell 3128 to 
~ 10 for A2163 with a median of ~ 2.5. The ratio of the 
number of AGN to the number of confirmed members pro- 
vides a strong upper limit to the AGN fraction in each cluster. 
If we use this larger value for the clusters with the smallest 
AGN fraction, rather than the value listed in Table [5] we still 
find the AGN fraction in these clusters is inconsistent with the 
values derived for the lower-redshift clusters with small sys- 
tematic errors. For example, the well-studied, high-redshift 
cluster AC114 has an AGN fraction of 5.2% based on con- 
firmed members alone and this upper limit is still inconsis- 
tent with the AGN fraction in clusters such as Abell 3125 and 
Abell 3128 that have small completeness corrections. The 
variation in the AGN fraction from cluster to cluster also does 
not seem to be due to radial dependence in the AGN fraction, 
as Figure [7] indicates that there is not a difference between 
the distributions of the Lx > 10 41 erg s" 1 population and other 
cluster members. 
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5.3. Correlations with Cluster Properties 

As there is evidence for variation in the AGN fraction from 
cluster to cluster, we compare the AGN fraction from each 
cluster with four properties that may influence the AGN frac- 
tion: cluster redshift; cluster velocity dispersion; amount 
of cluster substructure; and the fraction of Butcher-Oemler 
galaxies in the cluster. Simple arguments about the nature 
of AGN fueling and known relations between AGN and their 
host galaxies in the field motivates our choice of these prop- 
erties. First, a dependence with redshift is expected because 
of the pronounced decrease in the s pace densit y of AGN from 
high redshift to the present (e.g., Osmer 2004). In addition, 
the massive, early-type galaxies in clusters should harbor su- 
permassive black holes that most likely grew by accretion as 
luminous AGN at earlier times. Statistical studies of X-ray 
point source overdensities toward clusters do show evidence 
of an increase in the AGN fraction in high-redshift clusters 
(Dowsett 2005), although these AGN have not been spectro- 
scopically confirmed to be cluster members. Nevertheless, 
these arguments suggest that the AGN population in clusters 
of galaxies should be higher at higher redshift, although they 
do not necessarily require pronounced evolution over the red- 
shift range of this sample. In particular, most of the AGN in 
these clusters are low luminosity (Lx ~ 10 41-42 erg s -1 ) and 
current cosmic downsizing models indicate that there is rela- 
tively little evolu t ion at these lum inosities at low redshift (e.g., 
lUeda et al.ll200l lHasinger et al.ll2005l) . In fact the upper-left 
panel of Figure [10] actually shows that the AGN fraction in- 
creases at lower redshift for the Lx > 10 41 erg s" 1 sample, 
contrary to these expectations. However, there are substantial 
correlations between redshift and several other properties of 
these clusters and therefore this apparent correlation may be a 
selection effect. The most important of the potential system- 
atic effects produced by these correlations are due to the two 
lowest-redshift clusters, the merging pair Abell 3125/3128. 
These two clusters possess the greatest substructure and also 
have the lowest velocity dispersions. The cluster sample for 
this survey was simply chosen from available Chandra data 
in the southern hemisphere to obtain a measurem ent of the 
AGN fraction in clusters (see iMartini et al.ll2006T) and these 
potential parameter correlations were not foreseen. 

Perhaps the most important correlation between two cluster 
properties in our sample is that between redshift and veloc- 
ity dispersion. If mergers play an important role in triggering 
low-luminosity AGN, we expect that clusters of galaxies with 
lower velocity dispersions will have a higher AGN fraction 
because the merger rate will increase. The second column of 
Figure [10] shows that some of the lower velocity dispersion 
clusters do have higher AGN fractions, although this trend 
is largely driven by Abell 3125/3128. Recent evidence for 
an anticorrelation between AGN fraction and velocity d isper- 
sion comes from the work of Popesso & Biviano (2006), who 
studied the AGN fraction as a function of velocity dispersion 
with SDSS spectroscopic data. These authors found that the 
average AGN fraction increased in environments with lower 
velocity dispersion. 

Finally, the last two columns of Figure [10] show the AGN 
fraction as a function of cluster substructure P(<5) and the frac- 
tion of Butcher-Oemler galaxies in the cluster. The AGN 
fraction may depend on cluster substructure if cluster-cluster 
mergers produce a substantial increase in the galaxy inter- 
action rate. In the above discussion of how the merger rate 
is lower in higher velocity-dispersion environments, we ne- 



glected the potential role of fast fly-by interactions on fueling 
AGN. These fast interactions have been studied as a mech- 
anism for driving g alaxy evolution in clusters ( Moo re et al.l 
Il996t iGned in 2003) and may also provide a mechanism for 
fueling AGN dLake l et al.lll998l) . While it is not clear if the 
rate of fly-by interactions will increase in clusters with sub- 
stantial substructure, our measurements do indicate that the 
two merging clusters have the highest AGN fraction. Surpris- 
ingly, these clusters also have among the lowest fraction of 
Butcher-Oemler galaxies and this produces an apparent trend 
that clusters with fewer starforming galaxies have a larger 
AGN fraction. While again this is based on a small sam- 
ple that consequently does not merit too much interpretation, 
this is contrary to our expectation that clusters with more star 
formation and more available cold gas would be those with 
higher AGN fractions. 

6. SUMMARY 

We have investigated the distribution of AGN in clusters of 
galaxies in order to study AGN and galaxy evolution in rich 
environments. Specific questions that motivated our study in- 
clude the following: What mechanisms fuel accretion onto 
supermassive black holes? Is the evolution of AGN in clus- 
ters different from AGN evolution in the field? We have taken 
a step toward answering these questions with a study of the 
relative distributions of AGN and inactive galaxies in clusters 
and how the cluster AGN fraction varies as a function of sev- 
eral cluster properties. Our main results are as follows: 

1. The most luminous AGN (Lx > 10 42 erg s" 1 ) are more 
centrally concentrated than cluster galaxies with similar ab- 
solute magnitude Mr < -20 mag, while less luminous AGN 
{Lx > 10 42 erg s" 1 ) have a similar distribution to other cluster 
members. This greater central concentration of the most lumi- 
nous members is contrary to our expectation that most AGN 
will be triggered in the outskirts of clusters. While this survey 
does not extend to large enough radius to sample the cluster 
outskirts, our data do demonstrate that the Lx > 10 42 erg s" 1 
AGN are more centrally concentrated in the cluster than inac- 
tive galaxies of the same absolute magnitude and these AGN 
are not dominated by galaxies that have not recently fallen 
into the cluster. We would see newly infalling galaxies in 
projection in the kinematic distribution; however, as the clus- 
ter AGN have a similar kinematic distribution to other cluster 
galaxies, this suggests that they are not significantly more or 
less likely to be on a radial orbit than a typical cluster galaxy. 
Cluster AGN and inactive cluster galaxies also similarly trace 
cluster substructure. Future observations that encompass a 
larger fraction of the virial radial will help to physically sep- 
arate AGN in the cores of clusters from any potential AGN 
population in the cluster outskirts. 

2. The AGN fraction in clusters varies between clusters 
to a greater extent than can be explained by Poisson statis- 
tics alone. Specifically, the estimated AGN fraction /a (Mr < 
-2Q;Lx > 10 41 ) in cluster galaxies ranges from 14% to 0.6% 
with a mean fraction of 5% per cluster. While there are some 
systematic uncertainties due to our estimates of the cluster 
completeness, the statistical and systematic uncertainties for 
several clusters remain inconsistent with one another with 
greater than 90% confidence. We examined if the variation in 
AGN fraction from cluster to cluster correlates with any clus- 
ter property and found weak evidence that the AGN fraction is 
higher at lower redshift, in lower velocity-dispersion clusters, 
in clusters with substantial substructure, and in clusters with 
a smaller fraction of Butcher-Oemler galaxies. Unfortunately, 
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FIG. 10. — AGN fraction as a function of redshift, cluster velocity dispersion, amount of cluster substructure, and estimated fraction of Butcher-Oemler 
galaxies. The AGN fraction for each cluster is listed in columns 9 and 10 of Tableland the error bars are Poisson 90% confidence limits. The top row contains 
AGN with Lx > 10 41 erg s _1 and the bottom row contains AGN with Lx > 10 42 erg s _1 . The horizontal, dashed line corresponds to the average AGN fraction of 
the eight clusters for each luminosity threshold. The cluster properties are provided in Table[T] 

our sample of clusters is relatively small and has some sub- 
stantial correlations between these properties. In particular, 
the two lowest-redshift clusters are also those with the most 
substructure and the lowest velocity dispersions (Abell 3125 
and Abell 3128). Therefore we cannot conclude which of 
these properties, or what combination of these properties, sets 
the AGN fraction in clusters. A larger sample of clusters that 
more evenly fills this parameter space is needed to determine 
what parameters set the AGN fraction in clusters of galaxies. 
These data could then be used to test models of AGN fueling, 
for example via the functional form of the scaling between 
AGN fraction and velocity dispersion. 

Together these two results provide evidence that the prop- 
erties of cluster AGN, like the properties of cluster galaxies, 
vary both within clusters and from cluster to cluster. An 
equally interesting question to answer is if the AGN frac- 
tion varies substantially between clusters and lower-density 
groups and the field. While evidence from spectroscopic sur- 
veys unequivocally answers yes to this question and demon- 
strates that AGN are less commonly found in clusters than the 
field, radio and now our X-ray observations indicate that the 
dependence of AGN fraction on environment may be much 
less stark when the AGN are selected at these other wave- 
lengths. While differences in the typical spectral properties of 
AGN as a function of environment does not affect the present 
comparison between and within clusters, they will complicate 
comparisons of the AGN fraction between the field, groups, 
and clusters of galaxies. A possible resolution is that the 
AGN fraction may be a function of morphological type and 
not just absolute magnitude. This i s known to be the case for 
spectroscopically-classified AGN dHo et al.ll9 97) in the sense 
that early-type galaxies have a higher incidence of AGN. As 



iLehmer et alj ([2007) measured a comparable X-ray selected 
AGN fraction in morphologically-classified early-type galax- 
ies, and these dominate in clusters but not groups, future mea- 
surements of the AGN fraction may need to consider mor- 
phological type, in addition to host luminosity and AGN lu- 
minosity. This may reflect the fact that morphological type 
sets the spheroid mass and consequently the black hole mass 
for a galaxy of a given total luminosity. 
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